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Abstract. We give a general multiplication-convolution identity for the multivariate 
and bivariate rank generating polynomial of a matroid. The bivariate rank generating 
polynomial is transformable to and from the Tutte polynomial by simple algebraic op- 
erations. Several identities, almost all already known in some form, are specialization of 
this identity. Combinatorial or probabilistic interpretations are given for the specialized 
identities. 
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1. A convolution- multiplication identity 

This paper originated in an attempt to understand the relationship between Identities 
7 and 8 in this paper. One way to do this is to find a hidden identity which contains 
both identities as special cases. Identity 1 and its specialization, Identity 2, provide 
candidates for such an identity. As the reader will see, the proofs of Identities 1 and 2, 
and perhaps the identities themselves, are trivial. Identities have a way of receding to a 
Zen state of triviality. In the other direction, we add meaning by giving combinatorial 
or probabilistic interpretations for many identities derived from the hidden identity. 

Let M be a rank-r matroid on the set E with rank function rk and L(M) be its lattice 
of closed sets or flats. If M has no loops (that is, elements of rank zero), its characteristic 
polynomial x(M; A) is the polynomial in the variable A defined by 

X(M;A)= J2 M6,A)A" k «, 

X-.XeL(M) 

where n is the Mobius function in the lattice L(M). If M has a loop, then x(M; A) 
is defined to be the zero polynomial. The size-corank polynomial SC(M; a;, A) is the 
polynomial in two variables, x and A defined by 



SC(M;a;,A)= ^ x w A Mk(A) . 



A: ACE 



The size-corank polynomial has a natural "polarized" multivariate generalization. Let 
x be a labeled multiset {x e : e € E} of variables or numbers, one for each element e in 
the ground set E of the matroid M. If A C E, define x A to be the product Ile-eeA Xe • 
This notation is used analogously, so that, for example, —x = {—x e : e G E} and 
{—x) A — Y[ e . £ £a(~ Xe )- I n addition, if {y e : e S E} is another multiset, then xy is the 



set {x e y e : e £ E} and (xy) A = Y[ e - e eA x eVe- The subset-corank polynomial SC(M; x, A) 
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is defined by 

SC(M;x,X) = x A \ Mk{A) . 

A: ACE 

The subset-corank polynomial specializes to the size-corank polynomial when we set 
all the variables x e to the same variable x. Variants of the subset-corank polynomial, 
usually defined for graphs, have been rediscovered many times. The original discovery 
is by Fortuin and Kasteleyn in statistical mechanics. An almost complete list can be 
made by merging the lists in the surveys of Farr [J] and Sokal |14| . For matroids, the 
subset-corank polynomial is due to R.T. Tugger and it is sometimes named after her (see 

Em). 

The nullity-corank polynomial, usually known as the rank generating polynomial, 
R(M; x, A) is defined by 

R(M;x,\)= zW- rk ^A r - rk ^. 

A: ACE 

The size-corank and nullity-corank polynomials are closely related: indeed, 

ar r SC(M; x, xX) = R(M; x, A) 

and 

SC(M;x,X) = x r R(M;x,X/x). 
Both polynomials specialize to the characteristic polynomial. Specifically, 

SC(M;-1,A) = X (M;A), R(M; -1, -A) = (-l) r x(M; A). (1) 

The Tutte polynomial is defined to be the polynomial K(M;x — 1,A— f). The three 
bivariate polynomials transform into each other by simple algebra. (Note, however, that 
we may need to divide; thus, problems may arise when we wish to set a variable to 0,) 
We shall choose the polynomial that gives the simplest or most general identities. Thus, 
despite the title, we shall work mostly with subset-corank or size-corank polynomials. 
We begin with a general multiplication-convolution identity. 

Identity 1. Let M be a matroid on the set E, x and y be multisets of variables labeled 
by E, and A and £ be variables. Then 

SC(M;xy,X0 = ]T A r " rk ( T > (-y) T SC(M\T; -g, A)SC(M/T; y, Q. 

T: TCE 

Proof. We use the fact that if B and A are subsets of E such that B C A, then the 
sum 

E (-i) |THB| 

T: BCTCA 

equals except in the case A = B, when it equals 1. Then 



3 



SC(M;xy,XO 



V X B X r -MB) y A^r-rk(A) 



B,A: B.ACE 



E (-D 



|T|-|B| 



T: BCTCA 



T: TCE 



E 



k(B) 



B: BCT 



E /n (Mk 



(T))-(rk(A)-rk(T)) 



4:TCACB 



£ A r - rk ( T )(-y) T SC(A/|T;^,A)SC(M/r;y,0. 



T: TCB 



□ 



Setting x e = x, we obtain the following specialization of Identity 1. 

Identity 2. Let x, y, A, £ be variables. Then 

SC(M;xy,A£) = ^ A" k ( T )(-y)l T ISC(M|T; -x, A)SC(M/T; y, 0- 

T: TCB 

Next we state, without proof, the analog of Identity 2 for the nullity-corank polyno- 
mial. 



Identity 3. 

R(M;xy,X0= £ A Mk < T > (-y)l T l- rk ( T )R(M|T; -a:, -A)R(M/T; y, £)• 



T-.TCE 



The algebraic operation dual to convolution is comultiplication. Thus, we can express 
the identities in this section in the language of coalgebras. However, we will wait until 
we have more than just a formal theory. 



***Proof of Identity 30 



2. Weighted sums of polynomials 



Identities 1, 2. and 3 specialize to several known identities. We begin with identities 
expressing the size-corank or subset-corank polynomial as a weighted sum of other poly- 
nomials. 



This proof is included so that one can easily checked the identity. As indicated in the text, it will 
be removed in the final version. 



Identity 4. 

SC(M;xy,® = ]T (x + 1) T y T SC(M/T; -y, fl. 

T: TCE 

Proof. Set y e = — y e , A = 1, x e = — x e , leaving £ unchanged, in Identity 1 to obtain 
SC(M;xy,0 - X! y T SC(M|T;x,l)SC(M/T;-y,e) 

T:TC£ 



^ (g + l) VSC(M/r;-|/,Q. 



T:TCB 

In the last step in the derivation, we used a multivariate version of the binomial identity: 
SC(M\T;x,l) = Yl II x b 

B-.BCT b-.beB 



= H (X e + 1) = {X + 1) ' 



e: eST 



□ 



Identity 4 specializes to more familiar identities. 
Identity 5. 

SC(M;x,0 = Yl (x + l) X x(M/X;Q, 

X-.XeL(M) 

SC(M;x,0 = Yl (x + l) lXl x(M/X;0, 

X-.XeL(M) 

where the sums range over all closed sets X in the lattice L{M). 

Proof. If T is not a closed set, then SC(M/T; -1, £) = and if T is a closed set, then 
SC(M/T; —1,0 = x(M/T;£). Thus the range of the sums can be restricted to closed 
sets. □ 



R(M;xy,\0 

X \B\-^(B) X r-MB) y \A\-rk(A) ^r-rk(A) 

B,A: B,ACE 



E ^ r 



r-rk(T)/_s|T|-rk(T) 



E (-d |t| - |b| 

]T (_X)|B|-*(B)(_A)' 



T: BCTCA 

rk(T)-rk(S) 



^(|A|-|T|)-( r k(A)-rk(T))^(r-rk(T))-(rk(A)-rk(T)) 

A: TCACE 

52 ^- rk(T) (-?;)l T l- rk ( T 'R(M|T; -A)R(M/T; „, £). 
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The bivariate form of Identity 5 is a fundamental identity of Tutte [15] . Tutte found it 
for graphs and Crapo [2] extended it to matroids. The multivariate form appeared in 
[llj . It has the following interpretation. The subset-corank polynomial encodes the rank 
function of the matroid M in the sense that the rank of a set A is r — d, where d is the 
degree of A in the monomial x A X d in SC(M;x, A). On the other hand, SC(M; x — 1 , 
encodes the collection of closed sets of M in the following way: a set T is closed in M if 
and only if the monomial ca; T occurs in SC(M; x — 1 ,0 with a nonzero coefficient c. (The 
nonzero coefficient is x(M/T;A).) Thus, Identity 5 gives an "algebraic transformation" 
of the rank description of a matroid to its closed set description. 



3. Random matroids 

We next consider identities which give interpretations of size-corank or subset-corank 
polynomials as expected values of enumerative invariants of a random submatroid. 

We begin with an identity which is an "order dual" of Identity 5. Setting x e = 1, 
£ = 1, y e = — y e , and leaving A unchanged in Identity 1, we obtain 

SC(M;-y,A) = £ A" k ( T VSC(M|T; zl, A)SC(M/T; -y, 1). 

T-.TCE 

Using the fact that SC(M/T; — y, 1) = 1 — y E ^ T , we obtain the following identity. 



Identity 6. 

SC(M;-y,X) = £ X r -^ X (M\T; X)y T (l - y) ^ T . 

T-.TCE 

The next identity, obtained by setting x e = — x e , y e = —y e , £ = lj leaving A unchanged, 
is a generalization of Identity 6. 

Identity 7. 

SC(M;xy,X)= ^ A r - rk(T) SC(M|T; x, X)y T {l - yf XT . 

T-.TCE 

Bivariate versions of Identities 6 and 7, stated for graphs and given interpretations in 
terms of random graphs, are known. See Welsh [17] . Specifically, Identity 6 is related to 
an identity in Vertigan's Oxford thesis [16] and Welsh has given a proof using random 
subgraphs in [17] . A version of Identity 7 was found by Grimmett [7] (see also [17]). 

The interpretations by random subgraphs (given in [17) ) generalize easily to inter- 
pretations by random submatroids. Let M be a matroid on the set E. We generate a 
random subset and hence, a random submatroid, of M by deleting each element e in E 
independently and at random with probability l—p e . We need two somewhat artificial 
definitions, designed to make the interpretations work. Let N be a rank-s submatroid 
of the rank-r matroid M. Then the normalized characteristic polynomial x {N; A) to be 
■^ r ~ s x(N; A). Similarly, the normalized size-corank polynomial SC^(iV; x, A) is defined to 
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be A r S SC(N; x, A). Since the probability that the subset T is chosen is p T (l — p) T , it 
is clear that the expected value of the normalized characteristic (respectively, normalized 
size-corank polynomial) of a random submatroid of M equal SC(M; — p, A) (respectively, 
SC(M;px,X)). ~^ 

In analogy with random graphs, one can develop a theory of random sets of vectors 
or points in a finite ambient space. This was done in Kelly and Oxley [9]. There are 
two choices for the ambient space: the (affme) vector space GF(q) d or the projective 
space PG(d — 1, q) over GF(g). Here, q is a prime power and GF(q) is the finite field of 
order q. As a matroid, the projective space PG(d — l,q) is a simplification of GF(q) d . 
In particular, the two ambient spaces have the same lattice L(d,q) of flats. The two 
ambient spaces give the same theory of random sets, more or less. Since probability 
theorists usually work with random matrices and vectors, we shall work with GF(q) d . 
We define a random set V(d,q,p) of vectors to be a random submatroid of GF(q) d in 
which each element is chosen with the same probability p. 



Lemma. 



SG(GF(q) d ; x, A) = £ P (A - 1)(A - q) ■ ■ ■ (A - q d ^){x + If 

j=0 \J/ q 



f>C(GF(<z)V,A) 



d=0 



III 

.d=0 



1 + zq d 
\zq d 



5> 

d=0 



where 



[d\] q = (l-q)(l-q 2 )---(l-q d ) and 



[d!] 



m q [(d-j) 



Proof. The formula for SC(GF((7) d ; x, A) follows easily from known counting formulas 
in finite vector spaces and Identity 5. It is also a special case of formulas in Mphako [T3] , 
The generating function of SC(GF(g) rf ; x, A) factors into the product 

oo ^ oo ^ 

J> -1)(A -«).-. (A -g-WjL £(* + r/^ 

.d=0 1 ' iq ->- n 1 



d=0 



The formula now follows from the g-binomial theorem, in the version stated in [5] , Section 
1.3. □ 



Similar formulas exist for SC(PG(d — 1, q); x, A). Simply replace the exponent qi in x ql 

by q>- x +q j - 2 H \-q+l. 

We next give two typical results about expected values of matroid invariants of random 
vectors. We shall need to assume some knowledge of critical problems (see, for example, 

mm)- 



Theorem, (a) The expected number D(d,p,q, s) of s-tuples of linear functionals on 
GF(q) d distinguishing a random set V(d,p,q) equals 

fc=0 W 1 
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(b) The expected number sp(d,p,q) of subsets spanning GF(q) d in a random subset 
V(d,p,q) of vectors equals 



fc=0 



Wf) (-1)* ? G)(1+P)' 



Proof. Since the number of s-tuples distinguishing a matroid M represented as a mul- 
tiset of vectors in GF(q) d is the normalized characteristic polynomial x (M; q s ) evaluated 
at q s , 

D(d,p,q,s) = SC(GF(q) d ;-p,X). 

When A = q", the infinite product in the generating function telescopes into a finite 
product. Thus, 

oo ,i fs — 1 



^2D(d,p,q,s) 



d=0 



n( i+ ^) 



Ld=0 



£(i-p) 



d=0 



[dl] 



By an identity attributed to Euler or Cauchy (see, for example, [6], p. 254) 

s-1 



H(i-V)=E(- 1 )" 

d=0 fc=0 

The formula now follows from replacing the product with the sum and expanding. 

For part (b), observe that SC(M\T; 1,0) is the number of subsets A in T such that 
r — vk{A) = 0, that is, A spans M. Setting A = 0, y e = p e , and x e — 1 in Identity 7, we 
have 

sp(d,p,g)=SC(GF(g) d ;p,0). 
From the generating function, in the form given in eqn (3), it follows that 



J]sp(d,p, q) — 
d=o [a -^ 



,d=0 L '• 



Expanding the product yields the formula. This formula can also be obtained by Mobius 
inversion on the lattice L(d, q). □ 



Identities 4 and 5 can also be given probabilistic interpretations. If M is a matroid on 
the set E, we construct a random contraction H by choosing a random subset T by 
choosing each element e independently and at random with probability p e and letting 
H = Al/T. A typical result is that the expected value of the characteristic polynomial 
x{H; £) of a random contraction is 



4. The motivating identities 

We end by deriving the two identities which motivated this paper. We begin by setting 
x = —x, y = —1, and A = 1, leaving £ unchanged in Identity 3. Doing so, we obtain 
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the following convolutional identity, explicitly stated by Kook, Reiner, and Stanton in 
[5] and implicit in the combinatorial construction of Etienne and Las Vergnas in [3] : 

R(M;x,0= R(M|T;x,-l)R(M/T;-l,0. (4) 

T: TCE 

We can restrict the range of summation to flats for the same reason as in Identity 5. 
Because R(M/T; —1,0 is zero if M has isthmuses (or coloops), the range can be further 
restricted to cyclic flats, that is, flats with no isthmuses. 
Next recall that 

R(M; x, -1) = R(M ± ; -1, x) = (-l)' £; ' _ ' r x(A^ ± ; -x), 

where is the orthogonal dual of M. Thus, we obtain a version of the "KRSEV" 
identity which expresses the rank generating polynomial as a weighted-convolution of 
characteristic polynomials. 

Identity 8. 

(-l)H(M; J2 (-l) lXl x((M\X) ± ;x) X (M/X;0 > 

X: X£L°(M) 

where L°(M) is the lattice of cyclic flats of M. 

Next, we set x = 1, y = — 1, leaving A and £ unchanged in Identity 1, rederiving a 
multiplication identity in [12] . 

Identity 9. 

X (M;A0= £ ^ HX) x{M\X-\) X {M/X-e}. 

X: XeL(M) 

Kook, Reiner, and Stanton [8J have given an interpretation of Identity 8 for graphs 
using pairs of colorings and flows. One can easily adapt their interpretation using the 
critical problem. We will just give the result, referring the reader to [10] , Section 4.7, for 
the necessary background. Let M be a rank-r matroid on the set E and suppose that 
Q is an r x \E\ matrix over GF(s) representing M and R be an (l-El — r) x \E\ matrix 
over GF(£) representing the dual of M. If T C E, an (s, t)-duet (u, v) with support T is 
a pair of row vectors, such that u an |i?|-dimcnsional vector in the row space of Q with 
all coordinates in T equal to and all coordinates in E\T not equal to 0, and v is a 
T | -dimensional vector in the row space of the of R\T 1 the submatrix of R consisting of 
all the columns labeled by T, with all coordinates not equal to 0. Then 

R(M;t,») = (-l) r ^(_l)|su P port(T)| ; 

where the sum ranges over all (s, £)-duets (it, v). 

Interpretations of Identity 9 for graphs and representable matroids can be found in 
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